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We study a spontaneous collapse model for a two-level (spin) system, in which the Hamiltonian 
and the stochastic terms do not commute. The numerical solution of the equations of motions allows 
to give precise estimates on the regime at which the collapse of the state vector occurs, the reduction 
and delocalization times, and the reduction probabilities; it also allows to quantify the effect that 
an Hamiltonian which does not commute with the reducing terms has on the collapse mechanism. 
We also give a clear picture of the transition from the "microscopic" regime (when the noise terms 
are weak and the Hamiltonian prevents the state vector to collapse) to the "macroscopic" regime 
(when the noise terms are dominant and the collapse becomes effective for very long times). Finally, 
we clarify the distinction between decoherence and collapse. 
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I. INTRODUCTION 



Schrodinger equation in the Hilbert space 



-2. 



Stochastic Schrodinger equations find fruitful appli- 
cations both in the theory of open quantum systems 

allRM 11 and within collapse models |E 011 H 113, 
liaJiallllllflllTTj. In the first case, the stochastic 
terms mimic the effect of the environment on the open 
system^. In the second case, non-linear and stochastic 
terms are added to the standard Schrodinger equation 
to reproduce — under appropriate conditions — the col- 
lapse of the wavefunction^ ; in this way, it is possible to 
describe within one single dynamical equation both the 
quantum properties of microscopic systems and the clas- 
sical properties of macroscopic objects, and in particular 
the outcomes of measurement processes. 

In this paper, we analyze with numerical techniques 
the dynamics described by the following stochastic 
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^ The advantage of this kind of approach over the standard one 
based on the reduced density matrix of the open system is 
mainly computational: in some physical situations, it turns out 
to be computationally easier to solve numerically the stochastic 
Schrodinger equation and to perform the appropriate averages 
over different trajectories, rather than solving numerically the 
equation for the statistical operator. 

^ More specifically, it has been shown that a particular choice of 
the stochastic terms exists such that: 1) the quantum properties 
of microscopic systems are not altered in any appreciable way; ac- 
cordingly, these models reproduce all experimental results known 
at present; 2) at the macroscopic level, the reduction mechanism 
becomes rapidly effective so that macro-objects are always lo- 
calized in space. Note that in the literature, little attention (see, 
^■g-i 13 ) has been given to the behavior of the stochastic dy- 
namics when the Hamiltonian and the stochastic terms do not 
commute; here we will extensively discuss the consequences of 
this noncommutativity. 
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which is nonlinear and, as a consequence, nonunitary; 
anyway, it is easy to check that it preserves the norm of 
\ipt)- (cz) denotes the quantum average (V't|cz|V't): while 
ax^ (Tz are the usual spin operators. Wt is a standard 
Wiener process: ¥.[dWt] = and E [dWf ] = dt, and 
the (positive) parameter 7 measures the strength of the 
coupling between the spin operator and Wt. Finally, 
the stochastic dynamics is defined on a probability space 
(fi, J^, P) with the filtration {Tt : t g [0,oo)} generated 
in a standard way by the Wiener process Wt. 

Eq. ^ describes the evolution of the spin- 
wavefunction of a 1/2 spin particle; the dynamics is 
non-trivial, due to the presence of two competing terms 
which do not commute: on the one side the Hamiltonian 
H = hujax , whose effect is to rotate the state vector along 
the X axis of the Bloch sphere; on the other side, the spin 
operator which tends to localize the state vector along 
one of its two eigenstates |-|-), |— ). 

We will begin our analysis by studying the equation for 
the statistical operator p{t) — describing the ensemble 
of states each of which evolves according to Eq. — 
which in this case can be solved exactly (section II); we 
will give a qualitative description of the evolution of \4't), 
based on the dynamics for p{t). We will next solve Eq. 

numerically, showing that a much richer quantity of 
information can be obtained: in particular, we will ana- 
lyze the regime in which the reduction of the state vector 
to one of the two eigenstates of occurs, and we will 
clarify the distinction between collapse and decoherence 
(section III); we will analyze the reduction probabilities 
(section IV), the delocalization mechanism (section V), 
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and the transition from the "quantum" (7 <C w) to the 
"classical" (7 :s> uj) regime (section VI). We will conclude 
our analysis with a comparison between the model here 
analyzed and space-collapse models (section VII). 



II. THE STATISTICAL OPERATOR 

The statistical operator p{t) — E [\ipt) {ipt\] obeys the 
equation of motion'^: 



— p(t) = -iuj [(Tx,Pit)] 



(2) 



which of course is of the Lindblad-type [ig. It can be 
solved analytically for any given initial condition; if we 
consider the matrix elements of p with respect to the 
eigenstates |+), |— ) of (Tz: 



-\P\- 
-\P\- 



\P\-) 
\P\-) 



X y- 
y — iz 1 



- iz 

' X 



we get the following equations for x, y and z: 



x{t) = -2ujz{t) 
m - -2-1 y{t) 
z{t) = -uj + 2ujx{t) 



2jz{t). 



(3) 



The complete solutions of equations Q are written in ap- 
pendix I; here we limit ourselves to discuss some general 
features. 

The equation for y(t) is, trivially, a decaying expo- 
nential; the two equations for x{t) and z{t) lead to the 
following second-order differential equation for z{t): 



z{t) + 2-fz{t) + Alu^ z{t) = 0, 



(4) 



which is the equation of a damped harmonic oscillator; 
accordingly, for any value of lo and 7 (different from zero) 
the off-diagonal elements of the density matrix decrease 
exponentially in time. This results would suggest that 
reductions occur for any value of 7 — even though, of 
course, we expect the localization mechanism to be less 
efficient for smaller values of 7 — but we will see in the 
following sections that the real situation is subtler. 

The time evolution of x{t) has the following interesting 
property: for any value of lo and 7 (different from zero) 
and for any given initial condition, one has 



x{t) 



(5) 



This result is interpreted in the following way: suppose 
that the initial (at time t = 0) state vector is, e.g.: 




(6) 



^ The dynamical evolution embodied in Eq. has been studied 
in detail in Ref. 



and assume that after a certain amount of time it has 
been stochastically reduced either to the state |-|-) or to 
the state |— ), with the correct quantum probabilities. 
As a consequence, the density matrix becomes almost 
diagonal, with the diagonal elements equal to 3/4 and 
1/4. 

Result states that, as time goes to infinity, the di- 
agonal elements must change and become asymptotically 
equal to 1/2: this implies that, however it has been ini- 
tially reduced, the state vector starts jumping between 
the two eigenstates, so that after a certain amount of 
time there is an equal probability of finding it either re- 
duced around |-|-) or around |— ), despite that, at the very 
beginning, there was a much higher probability to find it 
reduced around |-f). In the following sections we will give 
a quantitative analysis of this effect. 

Note that the above result does not come as a surprise, 
since a general theorem by Spohn [T^ assures that, in the 
finite-dimensional case, when an equation of the quan- 
tum dynamical semigroup type admits a steady solution, 
then any other solution converges to it for t —^ +00. 
In our case, equation |2l has i/ as the only steady so- 
lution, where / is the identity matrix: accordingly, the 
off-diagonal elements must necessarily vanish as time in- 
creases, while the diagonal elements must converge to 
1/2. 



III. REDUCTION TIMES 

The dynamics described by eq. has several interest- 
ing properties. The first quantity which we have analyzed 
is the reduction time, i.e. the time it takes for the state 
vector to reduce to one of the two eigenstates |+), |— ) 
of the spin operator ct^ . A reduction process must satisfy 
the following two requirements: 

1. The state vector \^pt) must get "close enough" to 
one of the eigenstates of the reducing operator, cr^ 
in our case; 

2. Once the state vector reaches such an eigenstate, 
it must remain close to it for a "sufficiently long" 
time interval. 

Let us make these two requirements more precise. When 
IV't) gets close to the eigenstate |+) (|— )), the square 
modulus K+IV't)!^ (respectively, |(— IV'OP) approaches 
the value 1: accordingly, the first requirement amounts 
to saying that either K+IV'OP or |(— IV'OP is greater than 
1 — £, where e is a parameter depending on the physical 
situation under study. Following the second requirement, 
when — let us say — it happens that K+I^At)!^ > 1 — e, 
a reduction has occurred only if this condition holds for 
a time interval of length r; once more, the suitable value 
for T depends on the specific physical system. 

According to the above arguments, we define the re- 
duction time as follows: it is the smallest time tr such 
that, for t G [tr,tr + T], the square modulus K+IV't)!^ (or 



|(— iV't)! ) is greater than 1 — e. 
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FIG. 1: Reduction time tr as a function of 7; the simulations 
have been made for 7 = 5,10,20,40,60,80,100 s"\ while 
LO = 1 s~^. The picture shows the average values and the 
standard deviations over 100000 trajectories. Above each bar, 
we have indicated the fraction of the total trajectories which 
have been reduced within the time interval [0, 2tt] seconds. 



Fig. n shows the results of the numerical simulations 
for the reduction time as a function of 7, with the initial 
state vector'' as in eq. © and w = 1 s~^. We have 
taken e = 1/100, while r has been chosen according to 
the following criterion: when 7 = and jV'o) = |+) or 
|— ), the state vector rotates between the two eigenstates 
of CTz with frequency lo. Accordingly, even if no reducing 
terms are present, the dynamics brings periodically \ipt) 
close to, e.g., |+), in such a way that K+lfAt)!^ remains 
greater than 1— e for a time equal to: 7r/2— arcsin(l— e) ~ 
0.141 s. We have then chosen r to be equal to 10 x [7r/2 — 
arcsin(l — e)], to distinguish a reduction from the effects 
of the oscillatory motion induced by the Hamiltonian H. 

The simulation shows that, as expected, the reduction 
time decreases for increasing values of 7. When 7 is suf- 
ficiently larger than w (7 > 60 s""'), the reduction time 
becomes a well-defined quantity: all trajectories are re- 
duced, and the variance associated to tr decreases as 7 
increases^. On the other hand, when 7 is of the same 
order of — or smaller than — w, only a few trajectories 
(or none) are localized and the variance associated to the 
reduction times, increases: in this regime, the reduction 
process progressively disappears. This result is particu- 
larly important in relation to the time evolution of the 




* In all simulations, wc have taken JBJ as the state vector at time 
t = 0. 

^ The mathematical proof that, for ■y ^ lu, the state vector is 
reduced to one of the eigenstates of the "reducing oper ator " (ctz , 
in our case), was first given in ref. 0. See also ref. |lfil for a 
martingale proof. 



50 60 
t(s) 



FIG. 2: a) Time evolution of the real (dotted line) and imag- 
inary (continuous line) parts of {+\^t){'4^t\—) , for a typical 
trajectory followed by the state vector. The parameter 7 has 
been taken equal to 0.05 s~^, while as usual lj — 1 . b) 
and c) show the average values over 5 and 100 different tra- 
jectories, respectively. Note the decoherence effect: as time 
increases, the imaginary parts pick up different random phases 
and interfere destructively, d) shows the (theoretical) average 
values, corresponding to the time evolution of the real (y) and 
imaginary (z) parts of the off-diagonal element (+|p(i)|— ) of 
the density matrix. 



off-diagonal elements of the density matrix. 

In the previous section, we have seen that, for any 
value of 7, the off-diagonal elements of the density ma- 
trix decrease exponentially in time. According to the 
numerical simulations of Fig. ^ such a damping has two 
different origins according to whether 7 is significantly 
greater than lu or not: in the first case, the damping is 
due to the localization process; in the second case, it is a 
consequence of a phase randomization (decoherence) of 
the different trajectories. 

This property of eq. is well represented in Figs. 
121 and 13 Fig. shows a typical trajectory for the 
real and imaginary parts of (-|-|-!/'t)(V't|~)i when 7 — 0.05 
s~^: the real part (dotted line) decreases in time, while 
the imaginary part (continuous line) oscillates and does 
not approach the value zero. Figs, ^jp and ^ show 
the average values over 5 and 100 trials, respectively; 
while the time evolution of the real part of (+|'0t)('0t|~) 
does not change in a significant way, the evolution of 
the imaginary part displays a typical decoherence effect: 
the trajectories pick up different random phases and in- 
terfere destructively. Decoherence is then responsible 
of the exponential damping of the off-diagonal element 
{+\p{t)\—) = E[{+\4>t){'4't\—)] of the statistical operator, 
as shown in Fig. [21. 
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FIG. 3: As in Fig. |5| but with 7 = 50 s^^. Only the imaginary 
parts can be seen, as the real ones go rapidly to zero. Note 
that the decoherence effect is almost absent: when 7 is suf- 
ficiently larger than uj, single trajectories share the property 
that both the real and the imaginary parts of (+|^t){V't|~) 
are (most of the time) close to zero. 



Fig. Oi, instead, shows the result of the numerical 
simulations for the (typical) evolution of the real and 
imaginary parts of (+\^pt) {''Pt\—) , when 7 = 50 s""'^: both 
parts are almost always close zero, as an effect of the 
reduction process. Figs. ^ and ^ show an average 
over 5 and 100 trials respectively, while Fig. shows 
the mean value {+\p{t)\-) = E[{+\^pt){A\-)]- in this 
regime, the damping of the off-diagonal elements of the 
density matrix is not due to a decoherence effect but to 
the reduction process, since (almost all) realizations of 
(+|'0t)(V't|~) are vanishingly small. 



IV. PROBABILITIES 

The second interesting quantity to analyze is the prob- 
ability that the state vector reduces toward either the 
eigenstate |-|-) or the eigenstate |— ). The results of 
the simulation are shown in table UJ When 7 is much 
larger than u, the probability that a reduction occurs ap- 
proaches the quantum probabilities of finding the value 
or — 1 in a measurement of the observable Uz, given 
the initial state vector ©: 



Prob. [reduction near 
Prob. [reduction near 



I (+1^0) I' = 3/4 
I (+1^0) I' = 1/4 



This is an expected result. In fact, by denoting at 
{+\ipt) and (3t = {—\ipt), one has from eq. JQl: 



7 (s-^) 


red. near 4-) 


red. near j — ) 


total red. 


100 


73879 (74%) 


26121 (26%) 


100000 (100%) 


80 


74076 (74%) 


25924 (26%) 


100000 (100%) 


60 


73748 (74%) 


26252 (26%) 


100000 (100%) 


40 


73142 (73%) 


26858 (27%) 


100000 (100%) 


20 


70411 (71%) 


29033 (29%) 


99444 (99%) 


10 


45862 (66%) 


23279 (34%) 


69141 (69%) 


5 


3432 (61%) 


2165 (39%) 


5597 (6%) 



TABLE I: Reduction probabilities as a function of 7. The last 
column shows the number of times (over 100000 trials) the 
state vector has been reduced within the time interval [0, 27r] 
seconds. The second and third columns show the numbers 
of reductions toward the eigenstate |+) (second column) and 
| — ) (third column). 



In the regime 7 3> w, one can 



and similarly for \f3t\ 
neglect w: 



jatl = |ao|'+4V7 / M{1 - \a,\')dW,. (8) 
Jo 

Since the integrand in the above equation is bounded be- 
tween and 1, one can conclude ||20i that |at|^ is a mar- 
tingale, which in particular implies that the expectation 



value E [ [at |^] does not change in time: 



E[|at|^ 



|ao| 



(9) 



d\at\ 



^2ulm[atl3t]dt + A^\at\''{l - \at\^)dWt (7) 



in other words, reductions (and we have seen that in this 
regime all trajectories get reduced) occur with the correct 
quantum probabilities. 

On the other hand, when 7 is not much larger than uj, 
the approximation a; = is no longer valid, and the mar- 
tingale property does not hold any more. As a matter 
of fact, the numerical simulation shows that, when 7 de- 
creases, the reduction probabilities approach the values 

Prob. [reduction near I +)] ~ 1/2 
Prob. [reduction near [—)] ~ 1/2; 

anyway, we have seen in the previous section that for 
decreasing values of 7, fewer and fewer trajectories are 
localized, so that the reduction mechanism progressively 
loses its statistical meaning and eventually disappears. 

As a final note, we observe that the term 
—2ujlm[atPt]dt in Eq. Q — which spoils the martin- 
gale property of [atl^ — appears because the standard 
Hamiltonian H ~ hujax and the reduction operator Gz 
do not commute. If they commuted, such a term would 
be zero and |at|^ would be a martingale for any value of 
7- 



V. DELOCALIZATION TIME 

Even when 7^0;, the fact that the Hamiltonian H 
does not commute with the reduction operator Uz has 
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an important effect, which has not always been stressed 
in the literature: once the state vector has reached one 
eigenstate, it does not fluctuate around it forever, as 
there is always the chance that it jumps to the other 
eigenstate of a^- When this occur, we speak of a delocal- 
ization. 

This behavior of Eq. can be understood by looking 
at the equations of motion for the two components at = 
(+|V't) and/?t = (-IV't): 

dat = —iuj/Stdt — 2'yat{l — \at\'^)'^dt 

+2at{l-\at\^)dWt, (10) 

dPt = -iu;atdt-2-f(3til-\f3t??dt 

-2ptil - m'')dWt; (11) 

for simplicity, suppose that the initial state vector \^q) 
is equal — for example — to |+) or, equivalently, that 
the state vector has been previously reduced to the state 
1+)- 

We first of all point out that, if the Hamiltonian H 
commuted with the localization operator cr^ , and in par- 
ticular if H were equal^ to cr^, there would be no delo- 
calization from the state 1+), as at = e~"^' and /3t = 
would be the stationary solutions of the corresponding 
equations of motion. Accordingly, the fact that, in our 
model, H and do not commute, is crucial for the de- 
localization process. 

Coming back to Eqs. H1U|) and (|ll|l . an easy way to 
explain the delocalization mechanism is to look at nu- 
merical simulations. The noise dWt has a Gaussian dis- 
tribution with zero mean and variance dt, so on the aver- 
age it takes with the same probability both positive and 
negative values. 

Anyway, it can happen that — during a certain time 
interval — dWt takes (almost) only positive (or negative) 
values: it is precisely in this conditions that a delocaliza- 
tion may occur, as it is exemplified in tablclTTl There, we 
have simulated a solution of the equations of motions, 
using the following values: uj = 1 s~^, 7 — 100 s^^, 
dt = 10~^ s, ao = 1 and /3q — 0. The first delocalization 
begins at the 43937th step: we see from the table that 
the noise takes only negative values, thus determining 
the collapse of the state vector to the other eigenstate. 

We now move to analyze the delocalization time. It 
has been defined in the following way: assume that the 
state vector has been reduced to the eigenstate |-|-), which 
amounts to saying that |(+|?At)P remained larger than 
1 — e for a time interval at least equal to r. If, starting 
from any time td subsequent to tr, the square modulus 
becomes smaller than 1 — £ for at least a time 
interval equal to r, we say that the state vector has been 
delocalized from the eigenstate |-|-) and we call td the 
delocalization time. 



Which amounts to saying that, in Eq. llUi . the terms —iuif^t 
should be replaced by —iujat and, in Eq. illl . the term —iuiat 
should be replaced with iuj/3t. 



Step 


l«l' = i(+lv^>l' 


dWt 


43937 


0.98595512 


-0.01479509 


43938 


0.97223624 


-0.02990961 


43939 


0.91184014 


-0.04996265 


43940 


0.87878894 


-0.01746342 


43941 


0.80806981 


-0.02214550 


43942 


0.75135818 


-0.01448057 


43943 


0.55903192 


-0.02797212 


43944 


0.32168310 


-0.02466131 


43945 


0.07469396 


-0.03096076 


43946 


0.06038961 


0.00334229 



TABLE II: Piece of numerical simulation of Eq. Q. The 
values of the parameters are: uj — 1 s~^, 7 = 100 s~^, dt = 
10~^ s, and \il>o) = |+)- The table shows that, when dWt 
takes negative values for a sufficient number of steps, the state 
vector is driven away from the state |+) and reduced to the 
state |— ). 




5 10 20 40 60 80 100 

7 

FIG. 4: The picture shows — as a function of 7 = 
5, 10,20,40,60,80, 100 s"^ — the fraction of all trajectories 
(which have suffered a localization) that have been delocalized 
within the time interval [i7.,27r] seconds. 



Fig. 0]shows the results of the numerical simulation for 
the delocalization times. We have considered all trajec- 
tories which have been localized within the time interval 
[0, 27r] seconds, and we have counted the number of them 
which have gone through a delocalization within the time 
interval [tr,27r], where tr is the reduction time for each 
single trajectory. As expected, for increasing values of 
7, the fraction of all trajectories which are delocalized 
before time t = 2-k seconds, decreases: once the state 
vector gets localized to any of the two eigenstates of , 
it fluctuates around it for a time interval which is longer, 
the greater the value of 7. 
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sufficiently large, the reductions occur with the correct 
quantum probabilities. 

The state vector then fluctuates around the eigenstate 
to which it has been reduced for a time interval that — 
statistically — is longer, the greater the values of 7; ac- 
cordingly, when Lo can be neglected, the system behaves 
like a macroscopic system, regarding the physical prop- 
erties of the spin along the z-axis: they always have a 
well-definite value. 



FIG. 5: The figures sfiow two typical trajectories followed 
by the state vector in the Bloch sphere, for 7 = 0.05 s~^ 
(Fig. a) and 7 — 5 s~^ (Fig- b). In the first case, the 
Hamiltonian H rotates the state vector along the x axis of 
the Bloch sphere, while the reducing terms slowly move the 
rotation plane towards the |-|-) and |— ) eigenstates. In the 
second case, the effect of the reduction terms is so strong that 
the state vector is rapidly driven towards the |-|-) eigenstate 
and then fluctuates around it, until it jumps to the other 
eigenstate |— ). The arrows mark the starting point of the 
trajectories. 



VII. COMPARISON WITH SPACE-COLLAPSE 
MODELS 

Collapse models have been proposed in order to solve 
the measurement problem of quantum mechanics, and 
in particular to provide an explanation to the fact that 
macroscopic objects always have a definite position in 
space: accordingly, in most proposals on spontaneous col- 
lapse the localization operators are functions of the posi- 
tion operators'^. In this section, we confront the dynamics 
described by Eq. Q with the properties of the collapse 
model proposed in ref. The stochastic Schrodinger 
equation of this model takes the form: 



VI. TRANSITION FROM THE 
"MICROSCOPIC" TO THE "MACROSCOPIC" 
REGIME 



We conclude the analysis of Eq. by giving a general 
picture of how the evolution oi\ipt) changes when moving 
from the "microscopic" (7 <C uj) regime to the "macro- 
scopic" (7 ^ Lo) one*". When 7 is vanishing small, the 
state vector simply rotates along the a;-axis of the Bloch 
sphere, guided by the standard Hamiltonian H = fkuax- 
As 7 increases — still remaining smaller than lu — \ipt) 
keeps rotating along the a;-axis, but the stochastic terms 
shift the rotation plane towards the two eigenstates of (j^ 
(see fig. (Pi)). 

When 7 becomes significantly greater than lu, the state 
vector does not rotate any more around the x-axis of 
the Bloch sphere: it gets rapidly reduced to one of the 
two eigenstate of Cz (see fig. (P)); moreover, if 7 is 



-j:Hdt + ^ j (fx (iV(x) - (A^t(x))) dVFt(x) 



Sx{N{y.) - {Nt{x))fdt 



\^t), (12) 



where (7Vf(x)) — (V't|A^(x)|V't) and H is the standard 
Hamiltonian for the system under study. The "average 
density number operator" A^(x) is defined as follows: 



a 



d-'ye 



(x-y) 



aHy,s)a{y,s), 
(13) 

a^y,s) and a(y, s) being, respectively, the creation and 
annihilation operators of a particle at point y with spin 
s; Wt{x.) is a family of independent Wiener processes: 
E[dm(x)] = 0, E[dVKt(x) dWt{y)] = S^^H^ - y) dt. 

Finally, the two parameters, a and 7, appearing in the 
model, have been taken equal to: 



In collapse models 1^, the parameter 7 is constant and much 
smaller than the inverse of the characteristic times of the free 
Hamiltonian, so that microscopic systems are weakly coupled to 
the stochastic terms; this is the reason why the 7 a; regime 
has been called "microscopic". On the other hand, it can be 
proved that macroscopic objects strongly interact with the 
noise terms: in fact, the effects of such terms on the constituents 
of the macro— system sum with each other, in such a way that the 
dynamics of the center of mass of the object is governed by an 
effective single— particle collapse— equation, with a 7 much greater 
than the inverse of the characteristic times of the free Hamilto- 
nian. It is precisely this strong interaction which is responsible 
for their classical behavior. Accordingly, the 7 2> tt) regime has 
been called "macroscopic" . 



10" 



7 



47ry 
a J 



10" 



(14) 

Eq. H12|l has the same structure of Eq. and, 
in all interesting physical situations, the Hamiltonian 
H does not commute with the localization operators 
a^(y, s) a(y, s): as a consequence the dynamics described 
by Eq. (|I2I) is, qualitatively, analogous to that described 



* In refs llEl Il6l , an energy-based collapse model has been formu- 
lated; in this case, the Hamiltonian commutes with the reduction 
operators, so the model is different from the one considered here. 
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by Eq. We comment on the implications of our 

analysis of the spin-model for the space-localization 
model given by Eq. H12(l . 

Reduction mechanism. In the case of a free object, 
the analog of the time uj~^ can be reasonably identified 
with the time T it takes for the spread Ax of the center of 
mass of the object to change appreciably, e.g. to double 
its initial value; in the case of a free particle of mass m 
(in one dimension), we have: 

A.it)^A.ito)f7^1^^t^ (15) 

so that Ax{t) = 2Ax{to) when t = T, with: 

^ _ VT2m[Ax{to)]^ ^ 
h 

In the above formula, we shall take Ax{to) ~ 10^^ cm, 
which is the typical localization length of the model (com- 
pare with eq. ifTH) '). 

According to the numerical values given to the param- 
eters a and 7, in the microscopic situation (when only 
one or a few constituents are present) we are in the sit- 
uation where the characteristic time T associated to the 
Hamiltonian is exceedingly smaller than the characteris- 
tic time associated to the reduction mechanism; for 
a particle like a proton, we have, e.g.: 

T ~ 10-^s < X-^ ~ 10"s. (17) 

In this situation, the standard quantum dynamics of mi- 
croscopic systems in not appreciably altered by the pres- 
ence of the stochastic terms. 

On the other hand, in the macroscopic case, as already 
mentionedjEin "amplification mechanism" occurs: it can 
be shown [q] that the frequency A associated to the re- 
duction of the center of mass of a rigid body is approx- 
imately given by the formula: A — NX, where N is the 
number of constituents of the object. For an object of 
mass equal to 1 gr (iV ^ A^Avogadro), we have: 

T - lO^^s > A-^ - lO^'^s. (18) 

Accordingly, the characteristic time T associated to a 
macroscopic object is many order of magnitudes greater 
than A~^. In this situation, it can be proved that^: 

a) The reduction time is extremely well defined and small 
(~ 10^^ s), much smaller than the perception time of a 
human being. All macroscopic objects are localized in a 
very fast way, and superpositions of different macroscopic 
states cannot persist. 

b) After the reduction as occurred, the wavefunction of 
the system is very well localized in space: for typical 
macroscopic objects its spread in position is ^ 10^^^ cm. 



Of course, the wavefunction cannot have a compact 
support, since the Hamiltonian (when it does not com- 
mute with the position operators) would immediately 
turn it into a wavefunction with "tails" going off to in- 
finity; anyway, the fraction of total wavefunction cor- 
responding to the tails is of the order of e~^° , which 
means that they can be neglected for all practical pur- 
poses. 

Probabilities. The reduction probabilities are almost 
exactly equal to the quantum probabilities and the out- 
comes of measurement processes follow the Born rule. Of 
course, the physical predictions of Eq. (|12l) are in princi- 
ple different from those of standard quantum mechanics: 
this opens the way to the possibility of testing collapse 
models, especially at the mesoscopic level where the time 
scales of the quantum dynamics and the time scales of the 
reduction mechanism become comparable. Anyway, such 
kind of tests are difficult to implement as it is difficult 
to control all sources of noise which mask the effects of 
spontaneous localizations. 

Delocalization process. We have seen in section IV 
that, when the spin- wavefunction has been localized into 
one of the two eigenstates of az, it will sooner or later 
be driven by the dynamics of Eq. into the other 
eigenstate and back. The dynamics of Eq. H12|l displays 
a similar behavior. 

After a macroscopic object has been localized in space 
(take its velocity equal to zero), there is a non zero prob- 
ability that it gets delocalized and, shortly after that, 
localized again in a different region of space. Anyway, 
the probability that such an event occurs is exceedingly 
small^" (practically equal to zero) , being proportional to 
the fraction of the wavefunction corresponding to its tails. 
Accordingly, when a macroscopic object gets localized, it 
stays localized practically forever. 

Note also that this feature of space collapse is common 
to standard quantum mechanics as well: any macroscopic 
object, being described by a wavefunction with noncom- 
pact support, has a (vanishing small) probability of go- 
ing through a quantum jump from a space region to a far 
away one. 
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® See ref. llTl and references therein for all calculations and details. 



And, in particular, many orders of magnitude longer than the 
age of the universe. 
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Appendix I: Complete solutions of the equation for 
the statistical operator 

The solutions of equations JSJ are given as follows: 
Over damped case: 7 > 2aj. 



x{t) 




K he 


yit) 


= 2/0 




z{t) 




h de 



(19) 



with: 



to 



CO 



b = 



72^(1 + 
2^2 



7/i 



7/i 



^^0 



2/i 
2/i 



2^0 



LU 

^ 2^ 
2^ 



7 / 



Critically damped case: 7 = 2(x). 



(20) 



x{t) 

y{t) 

z{t) 



[a + 6t]e-T* 
2/oe-27t 

[c + dt] e-^\ 



1/2 



(21) 



with: 



a 

b 

c 
d 



(2x0 - l)t^ - 7^0 
(2x0 - l)t^ - 7^0- 



(22) 



Under damped case: 7 < 2cj. 



x(t) 
z(0 



[acos(Ai) + &sin(At)]e + 1/2 
[ccos(At) + dsin(At)] e~T^ 



(23) 



with: 



b = 

c ~ 

d = 

A = 



r"""2A 




(2cjxo — iij — 1 zq) 



(24) 



IE[l<+l-iC,>tf 
0.8 




FIG. 6: The two dotted curves represent the (analytical) evo- 
lution (+|p(i)|+), with 7 = 100 (upper curve) and 7 = 5 
(lower curve). The continuous curves show the average 
value of K+l''/'t)P c>v6r 100000 trajectories of lipt), for the same 
values of 7. For the values of the time step, see the text. 



Appendix II: A note on the numerical simulation 

To solve eq. (QJ, we have used the simple Euler al- 
gorithm, which proved to be satisfactory for the pur- 
poses of the article. To test the goodness of the sim- 
ulation, we have compared^^ the analytical evolution 
of E [K-hlV't)!^] = with the average value of 

K+lV't)!^ over many sample trajectories, obtained by 
solving eq. ^ numerically (see Fig. 

While analyzing the data, we have seen that for large 
values of 7 — i.e. when the reduction mechanism be- 
comes extremely rapid — the numerical simulation was 
not as good as expected. For such a reason we have de- 
cided to decrease the time step during the first 0.1 s, 
obtaining a sensible improvement. 

We have made different tests on the choice of the time 
steps and number of trials and, taking also into account 
the time required by our CPU for elaborating the data, 
we have found that the best choice (according to the pre- 
viously mentioned criterion) was to take: 

Time step for t < 0.1 s; 10""^ 
Time step for t > 0.1 s: 10"^ 
Number of trials: 100000. 



Of course, such a comparison tests the convergence of t he nu- 
merical algorithm in the weak sense (i.e. in the mean) |21| . while 
a test for the convergence in the strong sense (for single tra- 
jectories) would have been better; but this kind of a test is not 
feasible. Note that Euler scheme is of order 1 for the convergence 
in the weak sense, while it is of order 1/2 for the convergence in 
the strong sense. 
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